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ABSTRACT 

We have previously shown that the hncwidth distribution in AGN can be accounted 
for by an axisymmctric broad emission hne region. In this paper we show that the 
hncwidth distribution changes with redshift and that these changes are dependent 
on Ho and go- We show that relatively small samples of AGN at high redshift with 
measured linewidth at half maximum can be used to distinguish between values of Hq 
and qQ. Furthermore larger low redshift samples can be used to distinguish between 
luminosity functions and hence different models of quasar evolution. 
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1 INTRODUCTION 

It is generally expected that the width, v, of the broad emis- 
sion lines in active galactic nuclei (AGN) will show some 
dependence on the luminosity, L, of the nucleus as it is rea- 
sonable to assume that in higher luminosity systems the 
emitting gas have higher velocities. It is known from obser- 
vations that in variable sources linewidths vary in this sense 
with luminosity, e.g. Civ in NGC4151 (Fahey, Michalitsianos 
& Kazanas 1991). However it is not clear that the same depe- 
dence of linewidth upon luminosity can b e extended across 



samples of AGN. Observ ations of samples (Wills et al. 1993 



^uchnarewicz et al. 1997 ) show some evidence for a v — L re- 
lation, with however a large scatter. Linewidths must there- 
fore be dependent also on another parameter. It is known 
from observations that several properties of AGN appear 
to be dependent upon the line of sight to the observer. In 
radio lou d systems the linewidt h shows a clear dependenc e 
upon R dWills fc Browne 1986|; [Wills fc Brotherton 1995|), 



the ratio of core to lobe radio power, which is an indica- 
tor of viewing angle. While it is not clear that the broad 
emission line regions (BLR) are the same in radio loud and 
radio quiet systems it is not unreasonable to expect that the 
linewidths in radio quiet systems will also be viewing angle 
dependent. While it is possible that other parameters play 
a role, the minimal hypothesis is that v — v{L,i,X) with 
the same functional form for all systems, but different for 
the lines at various wavelengths A (since these have differ- 
ent widths). The viewing angle, i, is the angle of the axis of 
the BLR to the line of sight, i.e. i — for systems viewed 
face on. Since it is difficult to measure the angle to the line 
of sight for individual objects with any confidence, this hy- 
pothesis has to be tested statistically. To do so requires an 



■ galaxies: Seyfert - quasars: emission lines - cosmology: 



assumption about the angular distribution and the lumi- 
nosity function for AGN. For the angular dependence it is 
natural to expand in multipole moments to as high an order 
as the data justifies. 

In a previous paper (R,udge & Raine 1998) we used a 



random distribution in angle (cut off beyond some angle 
i *), th e luminosity function from Boyle, Shanks & Peterson 
( [l988| ) and a dipolar dependence on angle to construct a line 



width distribution function which could be matched to ob- 
servations. In that paper systems were generally at low red- 
shift and so cosmological effects were not considered. These 
cosmological effects enter in two ways. First through the lu- 
minosity function, and second through the minimum observ- 
able luminosity in the sample at each redshift z. Compar- 
ing the linewidth distribution at different values of redshift 
therefore provides, in principle, a method to determine the 
cosmological parameters Ho and qo. The aim of this paper 
is to test the feasibility of the method in terms of the num- 
ber of linewidth measurements required to determine these 
parameters. In fact, since we do not have an observational 
sample to test, we have to construct one; this requires as- 
sumptions about the cosmology. We shall therefore consider 
the related question of the number of linewidth measure- 
ments required to distinguish between significantly different 
cosmologies, between Hq = 50, 75 and 100 and between 
qo — 0.0 and go ~ 0.5. 

The method depends on knowledge of the luminosity 
function which itself depends on the cosmological parame- 
ters. This raises several issues. First, whether the method is 
any better than using fits to the luminosity function alone 
to determine the cosmology. While in principle this would 
be possible, in practice the dependence is weak and it is 
not easy to distinguish between different cosmological mod- 
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els from the luminosity data alone (Boyle et al. 1994). The 



implication is that the dependence of the linewidth distri- 
bution on the cosmology that arises through the luminosity 
function is also weak. However, there is a strong dependence 
of the shape of the linewidth distribution on the cut-off lumi- 
nosity which is itself sensitive to the cosmology. This has the 
following consequences for the way the method is applied. 

In principle we are seeking a solution for the cosmolog- 
ical parameters from simultaneous fits to the AGN luminos- 
ity function and linewidth distribution. However, the weak 
influence of the cosmology on the luminosity function means 
that we can use an iterative approach. We assume values for 
Hq and go in the luminosity function <^{L,z). We calibrate 
the linewidth distribution at low redshift and then compute 
it at high redshift for this <1>. Fitting to the data (when it is 
available for sufficiently many systems) gives a new Ho and 
go, which we can use to find a new $. This procedure can be 
iterated to convergence. In fact, the final effect of the weak 
dependence of the luminosity function on the cosmology is 
that very little (if any) iteration is necessary. 

A potential problem is that the luminosity function de- 
pends on intrinsic evolution as well as the cosmology. How- 
ever, we find that this shows up in an identifiable manner in 
the linewidth distribution, so it can be allowed for. Equiv- 
alently, we can use the method to solve simultaneously for 
the cosmology and to constrain the AGN evolution. We con- 
clude that to distinguish between open and closed cosmolo- 
gies (even if the Hubble parameter is not taken as fixed 
from other observations) would require in the range 100 to 
500 line profiles at a resolution of 500 km s~^ at a redshift 
of order 2. 



2 DEVELOPMENT OF THE z DEPENDENT 
LINE WIDTH DISTRIBUTION 



2.1 Line width distribution 



As in our previous paper (Fludge & Raine 1998) we assume 
a model for the BLR in which the FWHM is dependent pri- 
marily upon ionising luminosity and viewing angle alone. 
While it is accepted that there are other parameters which 
may have some effect on linewidth we have shown that such 
a restriction still allows us to account for the linewidth dis- 
tribution. Furthermore, in an axisymmetric model of the 
BLR, it is not unreasonable to expect these to be the most 
important physical parameters affecting linewidth. Thus the 
FWHM, V, of a given broad emission line is taken as a func- 
tion of the ionising luminosity and the inclination the sys- 
tem. This function can be expanded in spherical harmonics 
with luminosity dependent coefficients. In Rudge & Raine 
( 1998 ) we showed that the distribution of linewidths in low- 
redshift systems could be reproduced if this function were 
taken to be axisymmetric and only the first two terms of 
this series were retained and the coefficients taken to have a 
common dependence on luminosity. The FWHM of a given 
emission line is then given by 



= (a + 6 sin i)L2i 



(1) 



where the constants a, b and a are chosen for each emission 
line, i is the angle to the line of sight of the axis of the BLR 
and 1/44 is the luminosity in units of 10** erg s~^. 



Since it is difficult to determine the line of sight an- 
gle in a given system, at least with any accuracy, we are 
lead to consider the linewidth distribution rather than the 
linewidths of individual objects. Assuming that the inclina- 
tion of AGN is random across the sky, the number of systems 
per unit velocity range at each v is given for objects at low 
redshift by 



^(") = I ^$(^44)dL44 



(2) 



where the luminosity function ^{Lm) gives the luminosity 
distribution. To extend this to high redshift systems we need 
to take account of the cosmological and intrinsic evolution 
in "1> and the redshift dependence of the range of luminosity 
covered in flux-limited surveys. 



2.2 Luminosity functions 

The shape of the predicted linewidth distribution will clearly 
have some dependence on the shape of the luminosity func- 
tion. It is known from observations that the quasar popula- 
tion evolves with redshift causing the luminosity function to 
change with redshift. Thus it is natural to expect that the 
shape of the predicted linewidth distribution curve must also 
change with redshift. There is still much debate over how 
best to model this evolution. Once a model for the evolu- 
tion has been chosen, the parameters used to fit the data 
are then also dependent on the choice of values of Ho and 
go. Thus we shall first assume a luminosity function and an 
evolution model, and test for differences when Ho and go are 
changed. We then test for differences when the luminosity 
function and evolution model are changed. 

The quasar population can evolve in a combination of 
two basic ways. First, by pure luminosity evolution, where 
only the luminosity of quasars evolves with redshift and the 
total number remains constant. Second, by pure density evo- 
lution, where the total number evolves with redshift. Cur- 
rent work appears to favour pure luminosity evolution but 
it should be noted (e.g. Boyle et al. 1994) that there is still a 
great deal of uncertainty and no obvious best choice model. 

There are several models for the luminosity function 
available including those of Boyle e t al. ( 1994), Pel (1995), 
Maccacaro et al. (1991), Hasinger (1998) and Boyle et al. 
( 1998| ) . To begin with we shall concentrate on that of Boyle 
et al. (1994) because of the large sample size, the detailed 



information given for different models of evolution and the 
provision of results for two extreme values of go (0.0 and 0.5) 
spanning the range of possibilities. Th e lum inosity functions 
of Hasinger (1998) and Boyle et al. ( 1998|) s hould provide 
an improvement on this, with Hasinger ( 1998 ) using higher 
quality ROSAT data and Boyle et al. ( 1998| ) using ASCA 
data. However neither of these analyses contain information 
on the effect of changing go. The number of objects in the 
ASCA sample is also small. More information should be- 
come available for both of these lum inosity functions in the 
near future. Maccacaro et al. (1991) use only the Einstein 
EMSS {Extended Medium Sensitivity Survey , Stocke et al. 



( 1991 )) objects as used by Boyle et al. (1994) and develop a 
similar broken power law model of the luminosity function. 
Perhaps the best alternative to the luminosity func tion o f 
Boyle et al. (1994), for our purposes, is that of Pei (1995). 
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Table 1. Luminosity function p aram eters used for 2 power-law 
models G and H of Boyle et al. (1994). 
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1.53 
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43.70 
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1.89 
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1.36 
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1.73 


1.45 


Errors 


±0.15 


±0.1 


±0.2 


±0.1 


±0.1 





This uses a combined sample of around 1200 sources and 
tests two different models for the luminosity function. Note 
however that this is an optical rather than an X-ray lumi- 
nosity function. While this may seem an advantage in the 
following work, which focus es on the H/3 optical emission 
line, Maccacaro et al. (1991) emphasize the incompleteness 
of optical surveys and question the methods used to correct 
for this. In developing the redshift dependent linewidth dis- 
tribution we shall use the lu minos ity function of Boyle et al. 
(1994) and use that of Pei (1995) to test for differences in 

the resulting distribution functions. 

The luminosity function given in Boyle et al. (1994) 



gives several possible fits to the observed data using different 
evolution models. We have chosen to use models G and H 



of Boyle et al. ( 1994 ) as they are the best pair of fits to the 
combined ROSAT and EMSS data using the same evolution 
model but different values of go • The models are given in the 
form of a double power law 



$(L) = 



T 12 



L< L* 
L> L* 



(3) 



The dependence on z is included by evolution of L* . Lu- 
minosity evolution rather than evolution of the space den- 
sity is assumed for the models G and H. Thus, the quasar 
luminosity evolves by scaling in the following way: 



L{z) ^L{z^Q){l + zf 

L(z) = L{Z = Q)ll + Zma: 



■2 ^ Zfnax 
Z ^ Zmax • 



(4) 



It is important to note that (|3|) is in the form given in 
Boyle et al. ( 1994 ) and is correct for the de-evolved luminos- 
ity function i.e. the evolution model is applied to give all sys- 
tems as if they were at z = 0. To introduce a z dependence 
into (||), i.e. to replace $(-£/) by $(!/, z), the 2-power law pa- 
rameterization needs an extra scaling factor of (1 + 2)^^* to 
ensure that $(L* , z) is constant over z. Note that $* is given 
in units of 10"^ Mpc"^ (10*'' ergs s"^)"^ and thus must be 
scaled by L%^{z = 0) before entry into (Bl). The parameters 
for the two models are given in table |l. Fig. ^ shows how 
this luminosity function evolves with z and how the evolu- 
tion changes with different value s of g o and Hq. Note for 
later reference that Boyle et al. (L993) suggested that the 



optical and X-ray luminosities are related in the following 

,,ra,r T ^ j-O.88±'0.O8 

way Lx oc L^^^ 



The luminosity function of Pei (1995) is given by 

/3 I lV^ 



where the evolution model defines 
= L.(l + 2) 



{l+a)g-(^-</2<)^ 



(5) 



(6) 



The parameters for the model fits to the luminosity function 
are reproduced in table ^. Fig. |^ shows similar information 
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Figure 1. Broken power law model of the luminosity function 
by Boyle et al. 1994. Both panels show < z < 3 with the plot 
moving to the right with increasing z. The upper panel is for 
go = 0, Hq = 50 and the lower panel is for go = 0-5i ^^^0 = 100. 



Table 2. Parar netcr s of the Exponential L^/* law luminosity 
function of Pei (|i99e|). Note that h = Ho/100. 



Parameter 



log(L*/LQ 
log(<I>,/Gpc~'') 



{h, qo, 
(0.5, 0.5, 



a) 

-0.5) 



{h, go, 
(0.5, 0.1, 



a) 

-1.0) 



1.02 ±0.10 
2.73 ±0.05 
0.93 ±0.03 
9.87 ±0.10 
5.33 ±0.17 



0.66 ±0.10 

2.77 ±0.05 
0.91 ±0.03 

9.78 ±0.10 
4.34 ±0.20 



to fig. ^, but in this case for the exponential fit rather than 
the broken power law fit of Boyle et al. (1994). 

We now have two models for the luminosity function at 
two different values of qo- These have both been produced 
under the assumption that Hq = 50kms~'^ Mpc~^. To use 
other values of Ho requires only an application of simple 
scaling factors to L and $(L). These are 



L(//o = 50) = L{Ho)(^^y 
^L,z,Ho) = $(L,^,Ho = 50) (:^)'. 



(7) 
(8) 



Equation (^) scales the luminosity in a universe with Ho 
50 to the value to be used in (^) where Ho = 50. Equation (| 
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where the luminosity distance d_L is given by 

1)] 90 > 



dL 



l)(Vl + 2go2 



90 







.(10) 



Thus for a known flux limit and redshift it is possible 
to calculate the luminosity distance, di, and hence the 
minimum observable luminosity at this reds hift. For the 
ROSAT observations used in Boyle et al. (|l994, F ~ 



X 10 ^^ergss 



cm~ . In the following work, we assume 
that the samples generated will be complete down to some 
flux limit. While any incompleteness, particularly at the low 
luminosity end, will have some effect on the shape of the 
linewidth distribution and consequently the application to 
cosmology, for a sample of a few hundred objects at high 
redshift it should be possible to generate a luminosity func- 
tion specific to that sample. It is likely that tailoring the 
luminosity function to the sample, whether complete or not, 
will improve the accuracy of the results. 



2.4 New form for the linewidth distribution 

The linewidth distribution is now redefined to include the 
redshift dependent luminosity function and integrated over 
z so that the number of systems. A*', at each FWHM, v, is 
given by 

N{V)^ f / * ^""^ ^$(L44,z)dL44d2. (11) 



Figure 2. Pel (1995) exponential model of the luminosity func- 
tion. Both panels show < z < 3 with the plot moving to the 
right with increasing z. The upper panel is for go = 0.1, Hq = 50 
and the lower panel is for go = 0.5, Hq = 100. 



scales the calculated value of ^{L, z) to the value required 
for the universe where Ho ^ 50. These conversions reflect the 
fact that a smaller value of Ho would imply a more extended 
universe, making luminosities greater but thinning out th e 
galaxies, making calculated densities less (Weedman 1986). 



2.3 Limits of integration 

The redshift dependence of the linewidth distribution arises 
in just 2 places (assuming that the BLR itself does not 
change explicitly with z). The first is in the luminosity func- 
tion as detailed above. The second is through the z depen- 
dence of the limits of integration in (^). This is partly ac- 
counted for in the intrinsic evolution used in the luminosity 
function e.g. equation (^. However the lower limit of inte- 
gration is also affected by the flux limit of the observations. 
At some value of z the faintest source we can see is brighter 
than the evolved lower limit; in the Boyle et al. model this 
occurs for r^tZ"""^^" > L^iu{l+zf. So at any value of z the 
lower luminosity limit is the maximum of the two minima 
(the observational limit and the evolved lower limit). 

The minimum observable luminosity L corresponding to 
a given fiux limit F in the standard Freidmann cosmology 
(e.g. Marshall et al. 1984) is obtained from 



3 TESTING 

We now develop a method to test the way in which the 
linewidth distribution is affected by changes in qo, Ho and 
the luminosity function. Any significant differences in the 
shape of the distribution curve will clearly become more ob- 
vious at higher redshift. Any test therefore requires a sample 
of FWHM from high redshift AGN, say 2.4 < z < 2.6. As 
yet there are no published samples with a sizable number 
of objects and FWHM meas urements at high redshift . For 
example, the RIXOS sample ( Puchnarewicz et al. 1997) has 
only 6 objects in the range 2 < z < 3. The aim of the cur- 
rent work therefore has to be to provide a case for carrying 
out such observations by showing that the size of sample re- 
quired is not prohibitively large. This problem is approached 
in the following way. 

First we produce a set of model parameters a, b and 
Q for a sample of objects at low redshift where differing 
cosmology and evolution models make little difference. This 
is done for both luminosity functions. The value of be- 
tween the two model distribution curves is minimized by 
choice of a and b. Then, for a chosen pair of values for qo 
and Ho and luminosity function, we use these parameters 
to produce a set of model distribution curves at high red- 
shift. Once observational data is available this would then 
be compared with the theoretical distribution curve, and qo 
and Ho would be altered in an iterative process until conver- 
gence was achieved. However, as this data is not available, 
at this stage we generate a further model curve for different 
values of the cosmological parameters. A small random sam- 
ple is then generated from the first distribution curve and a 
test carried out to find how large this sample needs to be 
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so that the second distribution curve differs from the first 
at the 95 per cent confidence level. 

We also consider whether the distribution curve is sensi- 
tive to the small changes in the luminosity function model on 
changing go- This is done by generating a distribution curve 
for say go = 0.0 with the luminosity function as for go = 0.5 
and comparing to the true case where go = 0.0 in the lu- 
minosity function. The fits of the luminosity function to the 
observed data currently show no significant improvement by 
changing go. We will show that the linewidth distribution is 
more sensitive to the value of go than is the luminosity func- 
tion and thus go can be determined with a smaller data set 
by this method. 

This test of the effect of changing Hq and go is carried 
out on the linewidth distribution produced using both lu- 
minosity functions. A comparison between the two sets of 
distributions can then be carried out to find which features 
of the distribution curve are most sensitive to changes in 
the luminosity function. We can also test how large a sam- 
ple would be needed to distinguish between the luminosity 
functions. 

Note that for the test the distribution curves are 
transformed into histograms with bins of width 500 km s~^ 
to refiect the expected resolution of observations. The ran- 
dom sample is generated from this by a Monte Carlo type 
method to produce a similar histogram. The resulting his- 
tograms are then normalized to contain the same number of 
total objects. Thus the test is not sensitive to the abso- 
lute value of N{v) only to the shape of the curve. The test is 
not then sensitive to the evolution assumed in the luminos- 
ity function when only a small range of z is used. However it 
will be seen that plots of the peak height of the distribution 
against z do show very clear differences when the evolution 
model is changed. Discarding evolution at this stage of the 
testing is not unreasonable as the absolute value of N{v) is 
difficult to find observationally with a high degree of accu- 
racy, since this requires a high level of completeness in the 
sample. It will be seen that the evolution model is testable 
separately by consideration of the predicted dependence of 
N{v) upon z. 



4 RESULTS 

4.1 Low redshift distribution 

For a low redshift s ample we have selected th ose objects in 
the RIXOS sample ( Puchnarewicz et al. 1997) with z < 0.5 
and measured FWHM H/3. This provides a sample of 54 ob- 
jects including some for which the FWHM has considerable 
uncertainty. Fig. ^ shows the linewidth distribution for these 
objects overlaid with the model curve for a = 2500 kms~^, 
b = 7000 km s"^ and a = 0.35 using the luminosity function 
of Boy le et al. (1994). A similar fit is produced with the Pei 
(1995) luminosity function using a — 1000, b = 10000 and 
a — 0.31. Note that i*, the angle beyond which the BLR is 



totally obscured to view, is fixed at 60° for all calculations. 



4.2 Evolution with z 

While the greatest differences in the distribution are natu- 
rally expected to be seen at the highest redshifts, in order to 




5000 10^ 

FWHM (kin s"^) 

Figure 3. H/3 linewidth distribution for the RIXOS objects with 
2 < 0.5. 




5000 10* 
FWHM [km 




5000 10* 

FWHM [km s"^] 

Figure 4. Linewidth distribution at < z < 0.5 (upper) and 
2A < z < 2.6 (lower) for go = 0.5 (solid) and go = 0.0 (dashed). 
All curves produced using the Boyle et al. 1994 luminosity func- 
tion. 



make the possibility of verification by observation realistic 
we have chosen to use values for z in the range 2.4 < z < 2.6. 
This is separated from z = by a large enough gap and 
covers a sufficient depth to provide a reasonable chance of 
finding enough sources. This range is also towards the upper 
end of that covered by the luminosity function. 

Fig. ^ shows how the distribution around z — 2.5 
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z vs Distribution width 



z vs Distribution width 





z vs V at Distribution maximum 



z vs V at Distribution maximum 





vs Maximum N(v) 



vs Maximum N(v) 





Figure 5. Change in the distribution of linewidths with z for 
go = 0.5 (sohd) and go = 0-0 (dashed) produced using the 
Boyle et al. 1994 luminosity function. The upper panel shows the 
FWHM of the distribution, the middle panel shows the linewidth 
at maximum N{v) and the lower panel shows the value of N{v) 
at the peak. 



changes between go = 0.0 and go = 0.5. At high redshift 
there is a clear difference between the two distributions with 
that for go = 0.0 having a peak at higher FWHM, but a 
smaller number of systems at that value due to a broader 
distribution. However at low redshift the distribution for 
gp = 0.0 has a narrower width and hence a higher peak. Fig. 

shows how three key features of the shape of the distribu- 
tion change with redshift. 

(i) The full width at half maximum of the model curve. 

(ii) The linewidth at the peak value 

(iii) The value of A^(i;) at this peak 

Fig. ^ shows the same information as fig. ^ except that 
the shape of the luminosity function is fixed as that for go = 
0.5 to test whether the changes in the linewidth distribution 
with go are caused primarily by the changes in the luminosity 
fun ction. Both figs. ^ and g are produced using the Boyle et 
al. ( L994|)_ luminosity function. 

Fig. |7| shows the change with redshift of the distribution 
width, peak position and peak height produced with the 
different luminosity functions, with go = 1.0 and Hq — 50. 
We have carried out a similar test on the distributions 



Figure 6. As for fig. |5|but with go fixed at 0.5 in the luminosity 
function. Comparison with fig. H shows that the dependence of 
the distribution on z is therefore not driven by changes in the 
luminosity function. 



produced at 2.4 < z < 2.6 to find how large a sample would 
be needed to rule out one luminosity function assuming that 
the other is correct. As the test is not symmetrical the 
results are as follows. If the Boyle luminosity function is 
correct then a sample of about 150 objects would be needed 
to rule out that by Pei. Conversely, if Pei has the correct 
function then about 250 objects would be needed to rule 
out the Boyle form. 

These first results show that changes in the value of Hq 
and go do significantly affect the shape of the linewidth dis- 
tribution. However we need the results of the test to show 
whether these differences can be observed with a reasonable 
size sample. For each pair of go. Ho values we have generated 
a random sample of linewidths and tested this against other 
pairs of values. Table ^ shows the minimum sample sizes to 
reject the second pair of values in each c ase at the 95 per 
cent confidence level using the Boyle et al. ( 1994 ) luminosity 
function. 



Encouragingly many of the values in the table are of 
order 100 which is not prohibitively large; hopefully samples 
of this size will exist in the near future. However this is still 
beyond the scope of published samples. Table ^ shows the 
same information as table |^ except we use the Pei luminosity 
function and go = 0.1 in place of go = 0.0. 
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z vs Distribution width 




z vs V at Distribution maximum 




vs Maximum N(v) 




Figure 7. Comparison of changes in linewidth distribution with 
z using the Pei (1995) luminosity function (solid) and the Boyle et 
al. (1994) luminosity function (dashed). Both distributions pro- 
duced using Hq = 50 and go = 0.5. Individual panels are as for 
fig-i 



Table 3. Minimum sample size of quasar linewidths for which 
the pair of qo, Ho values in the first column can be rejected when 
the sample is generated using the pair of values in the first row 
of the table. Distributions were generated using the Boyle et al. 
1994 luminosity function. 
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5 DISCUSSION 

These results show clearly that the li newid th distribution 
changes with redshift. Mcintosh et al. (1998) find some ob- 



servational evidence for an increase in FWHM H/3 between 
t heir s ample at 2 < z < 2.5 and that of Boroson and Green 
( [l992| ) at low redshift, z < 0.5. This dependence of FWHM 
upon z may be biased by selection of only higher luminosity 



Table 4. As for table ^ except in this case the distributions were 
generated using the Pei (1995) luminosity function. 
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objects in the high redshift sample. More importantly this 
work shows that the evolution depends significantly on the 
values of Ho and qo and also on the choice of luminosity 
function. The test shows that we only need samples of 
~ 100 systems at z = 2.5 to be able to distinguish between 
extreme values of the cosmological parameters for a given 
luminosity function. While it is encouraging that these sam- 
ple sizes are not prohibitively large, the amount of FWHM 
data required is still beyond that currently published. There 
is however hope that the required data will be available in 
the near future. Searches for clusters of qu asars at high red- 
shif t, suc h as those by Boyle et al. (1997) and Newman et 
al. 



(1997) 



may well provide the necessary information. 
Perhaps the most significant result is that the major 
dependence on go is not within the luminosity function, as 
can be inferred from fig. ^. The X-ray luminosity function 
of Boyle et al. ( 1994 ) uses a large sample (~ 500) of ob- 
jects but cannot distinguish successfully between values of 
go • Using a linewidth distribution model should enable us to 
distinguish between values of go with samples of around 100 
objects. Current observational efforts to produce large sam- 
ples of quasars (e.g. Boyle et al. 1997) will provide us with a 
good source of candidates to observe further to obtain high 
resolution spectra from which FWHM measurements can be 
made. 

Comparison of luminosity functions is perhaps easier 
than testing for go and Ho as initially only low redshift sam- 
ples are needed. These sample do need to be much larger 
(~ 500) but there is much more available data at low red- 
shift. It may be that with a sample of this size, it is not 
possible to get a good fit to the linewidth distribution when 
certain luminosity functions are used. The distribution when 
using the Pei function was noticeably narrower than that for 
the Boyle function. If this does not prove possible, the large 
difference in the number of systems at the distribution peak, 
across the full range of redshift, should be a simple test of 
the evolution models with a small, but complete, high red- 
shift sample. Note that this statistic is not currently used in 
the testing due to the requirement of complete samples. 
In fact the distribution curves are all normalized to have 
the same total number of objects for the testing. Thus it is 
purely the shape of the distribution that is tested and not 
the actual number of systems. 

It should also be noted that current problems with the 
luminosity function models include not only the accurate 
shape and evolu tion but also t he luminosity range. It has 
been suggested (Hasinger 1998) that the AGN luminosity 
function should extend to low luminosities and join the nor- 
mal galaxy luminosity function. Our results should be sen- 
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sitive to this. Both of the luminosity functions used here are 
modeUed with a specific low luminosity cut-of f. 

In this work, as in Rudge & Raine ( 199^ , we have as- 
sumed that the distribution of sini is uniform and that the 
cut off angle i« is constant. In reality this is probably not the 
case. It is reasonable to expect that i* increases with lumi- 
nosity, i.e. the opening angle of AGN is greater in higher lu- 
minosity systems. In Rudge & Ra ine 1999 (unpublished) we 
have found in the RIXOS sample ( Puchnarewicz et al. 1997) 
that this is in fact the case, while the distribution of sin i is 
uniform in each luminosity bin. As a result the distribution 
of sini will be biased toward face-on objects. While the ef- 
fects that this has on the cosmological predictions are not 
clear as yet, any such effect will be reduced by removing the 
lower luminosity objects from samples. This, to some extent, 
is done naturally at high redshift as we cannot observe the 
lower luminosity sources assumed to be present. However re- 
ducing the range of luminosity considered will probably also 
result in an increase in the sample size needed. 

With the exception of the recently developed super- 
novae observations, all methods for obtaining qo suffer from 
uncertain evolutionary effects which introduce a scatter 
comparable to the magnitude of the effect being measured. 
Our method is no exception. What appears to be surprising, 
however, is that not only can the data be used to model the 
source evolution, but that this can be done in many cases 
with a relatively modest number of observations. 

While our work has concentrated on modelling the 
spread of linewidths by an axisymmetric BLR other authors 
have suggested other parameters than viewing angle to ac- 
count for the scatter in the linewidth-luminosit y rel ation, 
e.g. V = v{L,ax) mod el of Wandel & Boiler (1998) and 



the model of Robinson (1995) which has linewidth changing 
with profile curvature. 

In fact our results do not depend on the assumption 
made here that the scatter in the v-L relation is a viewing 
angle effect: any analysis which attributes this scaater to a 
single additional parameter will give similar results. Thus, 
while this method of obtaining the cosmological parameters 
is unlikely ever to achieve the accuracy or robustness of the 
supernova method, we believe it to be worth pursuing as a 
viable supplementary approach. 



6 CONCLUSION 

In this paper we have shown that in an axisymmetric broad 
line region, the linewidth distribution exhibits a strong de- 
pendence on Ho and qo- More importantly we have shown 
that these dependencies should be testable in the near future 
with observed samples of ~ 100 objects at high redshift. We 
have also shown that the changes are not primarily a result 
of changes in the luminosity function due to uncertainties 
in the value of go- It should also be possible to distinguish 
between model luminosity functions using large, low redshift 
samples or small, but complete, high redshift samples. This 
provides the theoretical basis for a new method of determin- 
ing the values of qo and Ho using observations of AGN. This 
statistical method makes no attempt to find the distance and 
luminosity of individual objects. 
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